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We theoretically investigate the dynamical Franz-Keldysh effect in femtosecond time resolution, 
that is, the time-dependent modulation of a dielectric function at around the band gap under an 
irradiation of an intense laser field. We develop a pump-probe formalism in two distinct approaches: 
first-principles simulation based on real-time time-dependent density functional theory and analytic 
consideration of a simple two-band model. We find that, while time-average modulation can be 
reasonably described by the static Franz-Keldysh theory, a remarkable phase shift is found to appear 
between the dielectric response and the applied electric field. 


I. INTRODUCTION 


In last two decades, intense coherent light of different 
characteristics has become available owing to advances in 
laser sciences and technologies. Ultrashort laser pulses 
can be as short as a few tens of attosecond, forming a 
new field of attosecond science [ij . Intense laser pulses of 
mid-infrared (MIIU or THz frequencies have also become 
available recently Employing these extreme sources 

of coherent light, it is possible to investigate the optical 
response of materials in real time with a resolution much 
less than an optical cycle [U, 0-[l|. 

The dielectric function e(a;) is the most fundamental 
quantity characterizing the optical properties of matter. 
Modulation of the dielectric function e{uj) in the presence 
of electromagnetic fields has been a subject of investiga¬ 
tion for many years. The change under a static elec¬ 
tric field is known as the Franz-Keldysh effect (FKE) 

, and that under an alternatin g el ectric field is known 
as the dynamical FKE (DFKE) [ITI - E^ . An important 
parameter which distinguishes DFKE from the static 
FKE is the adiabaticity parameter 7 = Up/fl, where 
Up = is the ponderomotive energy, and Ul 

is the frequency of the field, fi is the reduced mass of the 
electron, and E is the electic field [l^. A multi-photon 
picture applies for 7 << 1 , and a static FKE picture is 
appropriate for 7 >> 1 . Laser pulses having 7 ~ 1 is 
an intriguing regime where novel and unobvious DFKE 
phenomena are expected. 

In previous investigations of DFKE, the main focus 
was on the modulation of the optical response averaged 
over times much longer than the o ptic al cycle, examining 
the time-averaged fine structure [l^ and shifts in exci¬ 
tation structures In the present paper, we exam¬ 

ine DFKE in time domain, with a resolution much less 
than the cycle of the applied optical field. The DFKE 
response in subfemtosecond time resolution is very rel¬ 
evant to ultrafast optical switching in the teraherz or 
even petahertz (10^^ hertz) domain [1, Q- A first exper¬ 
imental report on the DFKE with a femtosecond time 
resolution has recently been given by Novelli et al. Q 
for GaAs, employing an intense pump pulse of THz fre¬ 


quency. They observed an interesting time shift between 
the pump pulse and the modulation of dielectric func¬ 
tion, but the mechanism of the observed time profile was 
not understood. To uncover the physics of time-resolved 
DFKE, we develop a pump-probe formalism in two dif¬ 
ferent theoretical approaches: first-principles numerical 
simulations based on time-dependent density functional 
theory (TDDFT [2^ ) and analytic investigation for a 
two-band model. Combining two approaches, we can un¬ 
derstand not only the strength of the modulation but the 
phase with respect to the pump field as well. 

The organization of the present article is as follows. 
In section II, we present formalism and results of our 
first-principles calculations for the time-resolved DFKE. 
In section III, we develop an analytical formalism for 
the time-resolved DFKE employing a parabolic two-band 
model. In section IV, a summary will be presented. 


II. FIRST-PRINCIPLES PUMP-PROBE 
CALCULATION 

A. Formalism 


In real-time TDDFT, we describe electron dynamics in 
a unit cell of a crystalline solid under a spatially-uniform 
electric field E{t). The method has been applied for 
calculations of linear optical responses 26 | and nonlin¬ 
ear electronic excitations by intense laser pulses 
Treating the field by a vector potential 


A{t) = -c / dt'E{t'), 


( 1 ) 


the electron dynamics in the unit cell of solid is described 
by the following time-dependent Kohn-Sham (TDKS) 
equation [2^ : 






+Vionir) + Vniut) -I- I4c(ut)] 


( 2 ) 


where me is the electron mass, Vion is the electron-ion 
potential for which we use a norm-conserving pseudpo- 
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tential [H, [s^, and Vnifjt) and Vxc{r,t) are electron- 
electron Hartree and exchange-correlation potentials, re¬ 
spectively. For the exchange-correlation potential, we 
employ an adiabatic local density approximation, using 
the same functional form of the potential for both ground 
state and time evolution calculations [s^- Since the 
Kohn-Sham Hamiltonian in Eq. ([J) has the lattice peri¬ 
odicity at each time, we may introduce time-dependent 
Bloch wave function, In prac¬ 

tice, we calculate the time evolution of the Bloch wave 
functions. 

We calculate electron dynamics in diamond, using a 
cubic unit cell containing eight carbon atoms. The TDKS 
equation is solved in real time and real space. The real- 
space grids of 22^ is used for the unit cell, and 32^ grids 
for the fc-points. The Taylor expansion method is used 
for the time evolution with a time step of At = 0.02 
in atomic unit. The number of time steps is typically 
70,000. An important output of the calculation is the 
average electric current density as a function of time, 
J{t). It is given by 


the absence of the pump electric field, the frequency- 
dependent conductivity a’(w) is related to the Fourier 
transforms of the electric field and the induced current as 
ct(w) = f dte^‘^*Jp(t)/ f Ep{t), where we may use 
any time profile for the probe electric field, Ep{t). We 
employ this relation in the presence of the strong pump 
field. To introduce the time Tp to explore the response, 
a pulsed electric field is employed for Ep{t) whose enve¬ 
lope shows a maximum at t = Tp. Namely, we use the 
following relation to define the time-resolved, frequency- 
dependent conductivity a(Tp,uj)^ 


d{Tp,uj) 


J dte'‘^*Jp{t) 

f dte^‘^^Ep{t)' 


( 5 ) 


Thus, in this consideration, the probe electric field has a 
dual role, to specify the time Tp at which we explore the 
response of the medium and to distort the medium to 
examine responses of the system. We consider that this 
dual role of the probe pulse will be adopted to measure 
the transient dielectric functions in most experiments. 


J{t) = -^ [ dr'^Retjj* (p+ in + JNL{t), 

rUeV Jv i \ c J 

( 3 ) 

where F is a volume of the unit cell. JNL{t) is the current 
caused by non-locality of the pseudopotential. 

It should be mentioned that A{t) in the TDKS equa¬ 
tion (I2|) is the vector potential in the medium and not 
that of the incident pulse in the vacuum. As we dis¬ 
cussed in Ref. , the relation between the two depends 
on the macroscopic shape of the materials as well as the 
direction of the polarization. Exchange-correlation ef¬ 
fects may also appear in the vector potential A^in time- 
dependent current density functional theory |4fll |. which 
we ignore for simplicity. 

To examine the DFKE, we carry out simulations solv¬ 
ing the TDKS equation (l2|) including both pump and 
probe electric fields in the vector potential A{t) [1^ . We 
assume that both pump and probe electric fields are lin¬ 
early polarized and are orthogonal to each other. We 
denote the pump electric field as Ep{t) and the probe 
electric field as Ep{t). The probe electric field is assumed 
to be weak enough to be treated by the linear response 
theory. We denote the electric current caused by the 
probe field as Jp{t), which is assumed to be parallel to 
the direction of the probe electric field. They are related 
by the time-domain conductivity as 

Jpit) = j dt'a(t,t')Ep{t'). (4) 

We note that the conductivity depends on both 

times t and t' rather than just time difference t — t' due 
to the presence of the pump pulse. 

We derive a frequency-dependent conductivity at time 
T from the conductivity We first note that, in 


B. Numerical results 


In the calculations below, we use the following elec¬ 
tric fields. The pump field is of the form Ep(t) = 
Eo,pfp{t) sinDt with the central angular frequency H set 
to fl = 0.4 eV; its direction is along the [001] axis. The 
field is turned on adiabatically described by the function 

/pW> 

fp{t) = I (0 < i < Tp) 

1 1 {t> Tp), 

( 6 ) 

where Tp is set to 10 fs. The probe field is of the form 

Ep{t) = Eo,p sm{ojpt) exp ^ , (7) 

oriented in the [100] direction. The average frequency ujp 
is set to 5.6 eV, which is equal to the calculated band gap 
energy of diamond in LDA. The field strength is set to 
Eq p = 2.7x 10“^ MV/cm, which is small enough to probe 
the linear response of the medium. The pulse duration rj 
is set to rj = 0.7 fs. With such a short duration, we may 
scan the spectrum of broad frequency region around the 
band gap from a single pump-probe calculation. 

The frequency-dependent conductivity is calculated 
using Eq. ([5]) with a slight modification. 


d{Tp,u}) 


f dte^‘^*g(t - Tp)Jp(t) 
f dte^‘^^Ep{t) 


( 8 ) 


where we introduced a filter function, g{t), [4l| to sup¬ 
press spurious oscillations which arise from a finite time 
period of the simulation. 
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FIG. 1. (a) Pump (blue-dashed line) and probe (red-solid 

line) electric fields are shown, (b) The electronic current 
(blue-dashed line) induced by the probe electric field (red- 
solid line), (c) The imaginary part of the dielectric function 
in the presence of the pump field, lm[e{oJ,Tp = 13/s)] (red- 
dashed line), and in the absence of the pump field, Im[e(ai)] 
(blue-dash-dot line). Black-solid line shows the difference, 
hxi[£{uj,Tp — 13/s)] — Im[£(a;)]. 


FIG. 2. Gontour plots (left) and their time averages (right) 
of AIm(£(a;, Tp)] under the MIR pump field of the intensity of 
(b,g) 5, (c,h) 7, (d,i) 10, (e,j) 20, and (f,k) 50 MV/cm. The 
time dependence of the pump electric field is shown in (a). 
In panels (a)-(f), the vertical solid (dashed) lines indicate the 
time of maximum (zero) of the pump electric field. 


Figure [T] shows an example of our pump-probe cal¬ 
culations. In (a), the electric fields of the pump field 
Ep{t) and the probe pulse Ep{t) are presented: The blue- 
dashed line shows the electric fields of the pump field 
Ep{t), and the red-solid line shows the electric field of 
the probe pulse Ep (t). The magnitude of the pump elec¬ 
tric field , Eq^p, is set to 20 MV/cm. The probe pulse is 
applied at a time when the magnitude of the pump field 
is maximum. The electric current induced by the probe 
pulse, Jp(t), is shown by the blue-dashed line in (b). The 
conductivity is calculated from the probe current using 
Eq. (IHl), and then converted to the dielectric function 
using a formula 

e(a;) = 1-I-47rf—(9) 
cu 

which is valid in the presence of the strong pump 
field. The imaginary part of the dielectric function, 
Im[e(w,Tp = 13/s)], is shown in (c). The red-dotted 
and blue dot-dashed lines present the dielectric func¬ 
tion with and without the pump field, respectively. A 
change of the imaginary part of the dielectric function, 


A Im[e(I7,Tp = I3/s)] =lm[e(0,Tp = I3/s)] — Im[e(u;)], 
is shown by the black-solid line. It indicates modulations 
of an exponential tail below and an oscillation above the 
band gap (w = 5.5 eV). 

Figure [5] shows changes of the imaginary part of the 
dielectric function caused by the strong pump field, 
AIm[e(w,Tp)] = Im[e(w,Tp)] — Im[e(w)]. In (a), time 
profile of the pump electric field is shown. In (b) - (f), 
changes of imaginary part of the dielectric function are 
shown in contour plots for four intensities. Horizontal 
axis is the time Tp and the vertical axis is the frequency 
w. In (g)-(k), modulations averaged over time are shown 
as a function of frequency. 

We first look at the case of strongest pump electric 
field, Eq^p = 50 MV/cm ,shown in (f) and (k). We find 
that an increase of the absorption below and a decrease 
above the band gap are seen when the magnitude of the 
pump electric field is close to the maximum. The modu¬ 
lation is small when the electric field is close to zero. This 
fact indicates the static FKE appears instantaneously fol¬ 
lowing the change of the pump field. We confirmed the 
modulation is well fitted by the static FKE formula if 
we assume the effective mass of /x = 0.25 to. With this 
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value of the effective mass, the adiabatic parameter 7 is 
29.5, much larger than unity, which is consistent with the 
appearance of the static FKE. 

As the field intensity decreases, we find changes in two 
aspects, one along the frequency axis and the other along 
the time axis. Above the band gap, we find an oscilla¬ 
tory behavior in the frequency direction. This is clearly 
seen in the time-average behaviors shown in (g)-(k). For 
example, at the electric field of 10 MV/cm, shown in (d) 
and (i), the modulation is negative between 5.6 eV (band 
gap) to 6.0 eV. It then becomes positive above 6.0 eV. 

Along the time axis, we can see a striking change in 
the phase between the modulation and the pump electric 
field. At the strongest electric field (f), the modulation 
and the pump electric field is in phase, as mentioned 
above. As the intensity of the pump field decreases, a 
phase shift forward in time is seen. The amount of the 
phase difference increases as the magnitude of the pump 
field decreases, as seen from (b) to (f). At the lowest in¬ 
tensity (b), the modulation signals appear at times when 
the pump field is close to zero. The parameter 7 for 
these intensities is (b) 0.29, (c) 0.58, (d) 1.19, and (e) 
4.76. Thus the phase change becomes appreciable when 
the 7 value is around and below unity. 


III. ANALYTIC CONSIDERATION 


A. General formula for the conductivity in the 
presence of a strong field 


the Houston function is expressed as 


=«nfep(t)(^ exp 


—i £ r 

I nkp{t') 


( 12 ) 


where kp{t) is defined by fcp(t) = k + eAp{t)/c. 

We then consider the solution of Eq. (fTOl) in the pres¬ 
ence of both pump and probe fields, A{t) = Ap{t)+Ap{t). 
We express the time-dependent Bloch function as 


“nk\ 


{r,t) = + (13) 


where the coefficients (t) are determined by the stan¬ 
dard procedure in the time-dependent perturbation the¬ 
ory. We have 

cLit) = f dt'PL{t')Ap{t') 

rUeC J_oo 

[ dt'kp{t')Ap{t'). (14) 

WeC J_^ 

Here we have introduced matrix elements of momentum 
operator. 




^^nn'kp{t) 

where 


exp 




,(15) 


We next develop an analytic consideration to under¬ 
stand behaviors of the DEKE signals seen in the numer¬ 
ical simulation. In the following developments, we con¬ 
sider a simplified description: electron dynamics in the 
presence of pump and probe fields is assumed to be de¬ 
scribed by a time-dependent Schrodinger equation for a 
single electron, 






-I- V{r) 


( 10 ) 


where V{r) is a time-independent, lattice periodic po¬ 
tential. We thus ignore the time-dependence of the 
Kohn-Sham Hamiltonian of Eq. except for the vec¬ 
tor potential. We express the solution of this equa¬ 
tion using time-dependent Bloch function as 

V'i(r,t) = 

We further assume that, in the presence of the pump 
field described by a vector potential Ap{t), the solution 
of Eq. (fTO is well approximated by the so-called Houston 
function [13, Using static Bloch orbitals u^^{f) and 
orbital energies which satisfy 


2 to{ 


-I- +V(r) 


"nk 


(r) = 




( 11 ) 


iPnn'k = f drKP^P^n'ki^ 
J 


(16) 


is the matrix element in the static basis. 

The electric current density averaged over the unit cell 
is given by 


J{t) 



Y. {Kk {p + ^^(^) + Imp) Wfe} ■ ( 17 ) 

nk 


This may be decomposed into pump and probe contribu¬ 
tions, J{t) = Jp{t) + Jp{t). The two components are 


Jpit) = 


e 

meV 


Y [ [M (^+ MP) , 

- Jv 

nk 

(18) 



^ lm{PYit)[PMt').Apit'))} (19) 


where Ue is the average density of valence electrons. 
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Equation (fT^ is a linear relation between the probe 
field and the induced electric current density. As we de¬ 
scribed in Eq. o, we introduce a conductivity function 
that relates the probe field and the induced 
electric current density, 





( 20 ) 


where a, /3 are the Cartesian indices of the current and 
the probe field, and and indicate a-components 
of probe electric current and probe electric field, respec¬ 
tively. The conductivity function for t > t' is given by 

) — ^e^a/3 


m'iV 


dt' 


’ 

n^n' ,k 


Im 


(Pa) 


nn'kp{t) 


iPd)n' 


n' nkp{t") 


X exp 



( 21 ) 


This is our general expression for the conductivity in the 
presence of a strong electric field, when the system is 
described by the time-dependent Schrbdinger equation 

. . . 

The polarization induced by the probe field is given as 
the integral of the current over time. We may introduce 
the linear susceptibility Xap{t,t') by 


pm= f dt'ut') = Y,f dt'xMt,t')E^p{t'). 

J —OO O 'J —OO 


_( 22 ) 

Therefore, the linear susceptibility and the conductivity 
are related by 



dt"aa0{t" ,t'). 


(23) 


B. Time-resolved, frequency-dependent 
conductivity 

As we described in Sec. Ill A[ we introduce a time- 
resolved, frequency-dependent conductivity by Eq. dSD 
using a probe electric field E^ (t) which has a sharp peak 
at t = Tp, 

^c,3(T„u;) = ^J—, (24) 

E^{uj) 

where and EP{uj) are Fourier transforms of J^{t) 

and EP (t), respectively. As a probe field, we first consider 
an impulsive probe field at t = Tp, E^(t) = k6{t — Tp), 
where k specifies the strength of the probe pulse. With 
this simple choice of the probe field, we have the following 
result for the time-resolved conductivity, 

dip{Tp,w) = J dse‘^‘^‘^aa/ 3 {Tp + s,Tp), (25) 


where the superscript I of indicates that the impul¬ 
sive probe field is used. As a more general probe pulse, 
we consider 


E^it) = fpit - Tp)e-“-(‘-^-) (26) 

where uip indicates an average frequency of the probe 
pulse and fp{t) is a real envelope function having a 
peak at t = 0. We assume that it is an even func¬ 
tion, fp{t) = fp{—t), so that the Fourier transform, 
fp{uj) = / dte^‘^*fp{t) is real. With this choice of the 
probe pulse and using Eq. (l24l) . we have 


- , Idsfpis)ai^iTp + s,uj)E(--^)^ 

j dsfp{s)ed---^)s ■ 


This expression indicates that the time-resolved conduc¬ 
tivity does not depend much on details of the shape of 
the probe pulse if we use a short enough probe pulse fp so 
that (Tp -I- s, uj) does not change much in the duration 
of the probe pulse. 

We note that the real part of the time-resolved con¬ 
ductivity Redaa(Tp,u;) is related to the energy transfer 
from the probe electric field to electrons, as the ordinary 
conductivity does. To show it, we consider the work done 
by the probe electric field to electrons which is given by 


W = 


dtEp{t)Jp{t). 


(28) 


This can be written as 


W = J2 [ dtEP{t)Jm 

a. 

= E^ / dioEP*iu^)J^{<^) 

^ 'J 



d^\El{u^)? 




du>\El^{L0)\^aaa{Tp,i0), 


(29) 


where we used the definition for the aapiTp,ui) given by 
Eq. (|27ll . Since the physical electric field, EP{t), and the 
physical induced current, JP{t), are real quantities, we 
have 


daaiTp,Uj) = aca{Tp,-Uj). (30) 


Using this relation, we have 



dw|TP(a;)pdaa(Tp,w) 



dw I TP (w) p Redaa (Tp, w). 


(31) 


This result clearly indicates that the energy transfer from 
the probe pulse to electrons is described by the real part 
of the time-resolved conductivity that we defined by Eq. 

(HZl). 
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C. Parabolic Two-Band Model 


We introduce a two-band model in Eq. (ED), consid¬ 
ering only two orbitals in the sum, occupied valence (v) 
and unoccupied conduction (c) bands. The excitation 
energy from the valence band to the conduction band is 
assumed to have a parabolic form, 


2/i 




(32) 


where Cg is the band gap energy and /i is the reduced 
mass of electron-hole pairs. 


1. Conductivity under a static electric field: static 
Franz-Keldysh effect 


For a static electric field E, the vector potential has a 
linear time-dependence, 

Ap{t) =Ap{T)-cE{t-T). (33) 


The conductivity a{t, t') of Eq. (I2T]1 is a function oit — t' 
and the conductivity a{Tp,uj) of Eq. (|25l) is independent 
of Tp. After straightforward calculations, we have 


rrieU) 




vck-eEs/2^P0) cvk+eEs/2' 




-(P/S) vck-\-eEs j2 cvk — eEs/2^ 


*{(efc+es)'S+“ 24 ^ 


(34) 


Below, we consider a real part of the diagonal element, 
Re(Taa(w)- We further assume that fc-dependence of the 
matrix elements may be ignored. Then, carrying 

out integration over s, we have 

Rea„„(u;) = -^|(p„),,|2 ^ j 

k 

X Ai |^(efc + eg-uj) ^ , (35) 


where Ai(x) is the Airy function. Further carrying out 
fc-integration, we have 


X 


Redaa (w) 

f eg-g; 

1 0 


(2/i)3/2e2 

2mluj 



\{Pa)vc\‘^VQ 



,(36) 


where 0 = (e^£'^/2/i)^/^. This is a well-known formula 
of static Franz-Keldysh effect 0. 


2. Conductivity under a periodic electric field: dynamical 
Franz-Keldysh effect 

We next consider a case of pump electric field which is 
periodic in time, Ap{t Tq) = Ap{t), where Tq is the 
period of the pump field and is related to the frequency 
^4 by Tq = 27r/f2. The conductivity has also the 

periodicity. 


aitA')=ait-Tn,t' -Tn). (37) 

We make a Fourier expansion of UapitA — s) which is 
periodic in t with the period To, 

00 

aa0{t,t-s)= e™*^‘cr^’^^(s), (38) 

n— — oo 

where (s) is defined by 

= - s). (39) 

The time-resolved frequency-dependent conductivity 
a^p{Tp,ijj) for an impulsive probe field defined by Eq. 
(l25l) may be expressed as 

^ipiTp, io) = Y (cc + nn), (40) 


where (cu) is the Fourier transform of (s). For a 
general probe field of Eq. (l26l) with the envelope function 
fp(t), the conductivity defined by Eq. (E7)) becomes 


da^(Tpj cc) — 


Jdsf{s)dipiTp + s,uj)C^^-^o^^ 
f dsf{s)ed‘^~‘^As 


„ /p(w - Wo) 


where /p(w) is the Fourier transform of fp{t). 

We now consider a specific form of the pump vector 
potential, Ap{t) = Aocosf2t, and calculate explicit form 
of the time-resolved conductivity assuming a two-band 
model. We calculate d^^p (w) in two steps. We first cal¬ 
culate the Fourier transform of crap{t, t — s) with respect 
to s which we denote as dap{t^uj). In the parabolic two- 
band model, it is given as follows. 


/ OO 

dse'^‘^^aa 0 {t,t - s) 

-00 


rUeU) 


-Tlp_6o, 




^00 _ 

I 




• ( 42 ) 
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The integral in the exponential is calculated as 

= (cfc + Cg + Up) s — 9i sin r 2 (t — s) + 0 i sin fit 
—02 sin 2fl(t — s) + 6*2 sin 2fl<, (43) 


where we introduced 


U„ = 


e^Al 




01 = 




ek ■ Aq 
ficH, ’ 


0 =^. 
^ Siic^Q 


(44) 


(45) 


(46) 


Using a relation involving Bessel function Jn{x), 

e“"in® = ^ J„(a)e“®, (47) 


where n runs for whole integers, we may express 
sin sin 2Q.t = E 


I'm 


= Y,MOi.02)^ 


AlQt 


(48) 


where we defined the generalized Bessen function by 

^ 2 ) = ^ d;_2m(^l)'^m(^2)- (49) 

m 

We may express the integral in Eq. (l42l) as 


f 


dsU'^’^ exp 


+01 sin fit — 9s sin2fl(t — s) + 022 flt] 


E LJ 

in — 


tJq( 0 i, 02 )J; 2 ( 0 i, 02 ) 


I 1 I 2 


(cfc + Cg + Up + lifl) 


(50) 


Equation (|42ll is now written as 




rrieOJ 


kUlo 


{Pa)vc {Pl3^c 

(,P0 )'uc {_Pa ) cv 




u! — (cfc eg + Up + lifl) 

0 + (cfc + Cg + Up + Zlfl) 


( 51 ) 


We next calculate (w) as follows. 


^ die ^^^*aap{t,uj) 


^a/3 ' 


• 9 -2 

'kle^O'BOnO 


rrieUJ 


mlwV 


E 

kl 


Jl{9l,92)Jl-n{9l, 02 ) 

— J;( 01 , 02 ) 0 Z+„( 0 i, 02 ) 


iPa)vc{pp)c 


OJ — (Cfe + Eg + Up + Zfl) 

{pp) vc{Pcx) cv 


^ + (cfe T Eg + Up + Ifl) 
We put the above expression into Eq. dm, 


.(52) 


frr \ /p(w + nfl — Wo) inQ,T„~{n), , 

^p{Tp,uj) = ^-e PcrW(a; + nfl) 

„ /p(w - Wo) 


te 


-TleSn, 


TTleW 


+ E: 


le 


fp{uj + nfl - Wo) 


klr 


\{u: + nn)V fp{uj-u;Q) 
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We carry out k integration in terms of Ck = k^/2p. and 
cos0fc. After integration over cos 9k, only even n terms 
contribute. We introduce the following quantity, 

6n(fc) = J d(cOS0fe)Ji(0i,02)Ji-n(6'l,02), (54) 

and we denote n = 2m below. For the diagonal term 
a = /?, we have 
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As the final result of this subsection, we obtain the 
following expression for the real part of the conductivity, 
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CoefRcients Cm{i^) are given by 
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where k and are related to Ck and by fc = 
are defined by 


= ±(a; + 2mVl) — (eg + Up + Ifl). (58) 


Coefficients S^i^) are given by 
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We note that the terms with m = ±1 are responsible for 
the phase shift seen in Fig. [2] 

The above expression of the time-resolved conductiv¬ 
ity depends on the choice of the envelope function fp(t) 
through fp{uj). Modulations with a long period are in¬ 
fluenced by fp{uj) at small w, typically w ~ ±217. This 
indicates that the conductivity does not depend much on 
the detail of the probe pulse, if we employ a probe pulse 
which is much shorter than the period T of the pump. In 
the TDDFT calculation shown in Fig. 2, we used a short 
enough probe pulse with /p(±2f2)//p(0) = 0.96 close to 
unity. 


D. Numerical results 

We numerically calculate the change of the dielectric 
function AIm[e(w,Tp)] using Eqs. (l56l) - (l5^ . In the 
calculations shown below, we simply put /(w) = 1. This 
is equivalent to using an impulsive field for the probe, 
Ep{t) oc S{t - Tp). 

Figure [3] (a)-(d) show the contour plot, where the hor¬ 
izontal axis is the phase defined by QTp. In Fig. [3] (e)- 
(h), modulations averaged over time is shown by red-solid 
lines. Blue-dashed lines show results of static FKE: the 
time-averaged modulation using static formula , Eq. (1351) 
where the parameter 9 is evaluated using time-dependent 
electric field, E{t). The parameters in the two-band 
model are set to reproduce responses in the TDDFT cal¬ 
culations: The frequency of pump field is set to D = 0.4 
eV , the reduced mass is set to ^ = 0.25me, and the 


(a) 

5 MV/cm 


3 


(b) 

1 0 MV/cm 


3 


(c) 

20 MV/cm 


3 


(d) 

50 MV/cm 


3 



FIG. 3. Contour plots and time averages of AIm[£(a;, Tp)] in 
the two-band model for the pump field intensities of (a,e)i?o = 
5, (b,f) 10, (c,g) 20, and (d,h) 50 MV/cm. In the left panels 
(a)-(d), the horizon axis is the phase defined by DTp. The 
vertical solid (dashed) lines indicate the position of maximum 
(minimum) of the electric field. In the right panels (e)-(h), 
red-solid lines is the time average of the modulation. Blue- 
dotted lines show result of static FKE. 


dipole matrix element is set to IPcuP = 0.928 in atomic 
unit. 

We first look at time-averaged modulations shown in 
panels (e)-(h). In all panels, we find appearances of ab¬ 
sorption below the band gap and decrease above the band 
gap. We can also find structures of at w = Cg ± D, most 
clearly in (b) and (f), which are cause by side band con¬ 
tributions with different I values in Eqs. eH) and (j59() . 
They show a similar behavior to Fig. [2] (g)-(k). We also 
find that the time-averaged modulation is quite close to 
the estimation by the static FKE, even at the smallest 
intensity where the 7 value is much less than unity. We 
repeat a similar calculation changing the frequency of the 
pump field while the intensity is fixed at 50 MV/cm, and 
have found that the static FKE describes reasonably the 
time-averaged modulations below the band gap. 

We next look at the time dependence of the modula¬ 
tion. At the strong pump electric field (d), we find a 
large modulation when the magnitude of the pump elec¬ 
tric field is large. We find a phase shift as the pump 
electric field decreases. The maximum of the modula¬ 
tion moves forward in time. From these observations, we 
can say that all the features seen in the first-principles 
TDDFT calculations shown in Fig. [2] are reproduced by 
the analytic formula of Eqs. ((56ll - ((59l) of the simple 
two-band model. 
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IV. SUMMARY 


We have developed a theoretical pump-probe formal¬ 
ism to investigate the dynamical Franz-Keldysh effect 
in time-domain, in femtosecond time resolution much 
shorter than the optical cycle of the applied pump field. 
Both numerical simulations based on real-time time- 
dependent density functional theory and an analytic ap¬ 
proach in the two-band model reveal the same behav¬ 
ior for the modulations of dielectric properties. We find 
the time-averaged behavior in the DFKE can be well de¬ 
scribed by the static FKE. The most remarkable feature 
of the DFKE is that there appears a phase shift between 


the modulation of the dielectric function and the applied 
pump field which becomes significant as the magnitude 
of the electric field decreases. 
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